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ABSTRACT 

LS I +61 303 and LS 5039 are exceptionally rare examples of HMXBs with MeV-TeV emission, 
making them two of only five known or proposed "7-ray binaries" . There has been disagreement 
within the literature over whether these systems are microquasars, with stellar winds accreting 
onto a compact object to produce high energy emission and relativistic jets, or whether their 
emission properties might be better explained by a relativistic pulsar wind colliding with the 
stellar wind. Here we present an attempt to detect radio pulsars in both systems with the Green 
Bank Telescope. The upper limits of flux density are between 4.1-14.5 iiJy, and wc discuss the 
null results of the search. Our spherically symmetric model of the wind of LS 5039 demonstrates 
that any pulsar emission will be strongly absorbed by the dense wind unless there is an evacuated 
region formed by a relativistic colliding wind shock. LS I +61 303 contains a rapidly rotating Be 
star whose wind is concentrated near the stellar equator. As long as the pulsar is not eclipsed 
by the circumstellar disk or viewed through thcj^densest wind regions, detecting pulsed emission 
may be possible during part of the orbit. 



Subject headings: pulsars: general - stars: individual(LS I +61 303, LS 5039) 



1. Introduction 

LS I +61 303 is a high mass X-ray binary 
(HMXB) that is also a confirmed source of very 
high energy 7-ray emission. The system consists of 
an optical star with spectral type BO Ve and an un- 
known compact companion in a highly eccentric, 



26.5 day orbit (|Aragona et al. 2009 


Grundstrom et 


al.||2007 


Casares et al.||2005a 


). While the system 



has a relatively low X-ray luminosity for a HMXB, 
LS I -1-61 303 is the 15th brightest 7-ray source in- 
cluded in the Fermi LAT 1-year Point Source Cat- 
alogue (Abdo et al. 2010). The Be disk interacts 



with the compact companion, producing orbital 
phase modulated emission across the electromag- 



netic spectrum: TeV (Albert et al. 


2006 20081, 


GeV 


(Abdo et al. 


2009a 1, X-ray (F 


'aredes et al. 


1997 


Leahy 2001 


Torres et al.||2010 


I, optical Ha 


(McSwain et al.|2010), and radio (Gregory & Tay- 


lor|1978 Taylor & Gregory|1982 Ray et al.|1997|. 


Taylor & Gregory 


( 1982 1 defined the arbitrary ref- 



erence for zero phase at JD 2,443,366.775 that re- 
mains the conventional definition for LS I +61 303. 
The binary orbital period (P = 26.4960 ± 0.0028 
d) has been determined from the periodic radio 
outbursts that peak near 0(TG) = 0.6 — 0.8 (Gre- 
gory|2002 ), and this period is confirmed by optical 
spectroscopy ('Arago na et al.|2(309 1 . Periastron oc- 
curs at (/)(TG) = 0.275 ( [Aragona et al.|[2009 ). 

LS 5039 is a HMXB that consists of a mas- 
sive, main sequence star (ON6.5 V((f)) spectral 
type; McSwain et al. 2004) and a compact com- 



panion in a tightly bound, eccentric orbit ( Casares 



et al. 


2005b 


2010 


Aragona et al. 


2009 


) . Similar 



to LS I +61 303, LS 5039 exhibits 7-ray (Abdo 



et al. 2009b Aharonian 2006) and X-ray (Bosch 



Ramon et al. 2005 2007 Takahashi et al. 2009) 



emission that is variable over the 3.9 day orbital 
period. There is also evidence for sub-orbital mod- 
ulation of the stellar radial velocities, possibly due 
to tidally excited nonradial pulsations of the O- 
type star ( Casares et al.|[20T0 ). 

The systems LS I +61 303 and LS 5039 recently 
became subjects of controversy regarding the ori- 
gin of their high energy emission and the nature 
of their compact companions. Many studies (e.g. 
Massi"2 004] [Paredes et al.||2002l [Bosch-Ramo"iret 
al. 2005, Pandey et al. 2007) have assumed the 



"microquasar" model - so named because they 
are observed with extended radio emission that re- 



semble small-scale versions of the relativistic jets 
of extragalactic quasars. [Yamaguchi fc Takahara] 
(pOTol) argue that the keV-TeV variability of LS 



5039 is due to to the injection of electrons into 
the relativistic jet in a low inclination binary orbit 
(which requires a black hole given the observed 
mass function; Aragona et al. |2009| [Casares "et| 



al. 20101. Romero et al. (2007) argue that LS I 



+61 303 is an accreting microquasar based upon 
Smoothed Particle Hydrodynamics models of the 
stellar wind and the variable accretion rate during 
the orbit. However, there is no observed evidence 
of an accretion disk in either system. Optical spec- 
troscopy may reveal Ha emission and enhanced 
continuum emission from an accretion disk, but 



neither are observed in LS 5039 (McSwain et al. 



2004). Emission from the Be star's circumstel- 
lar disk dominates the Ha line profile in LS I +61 
303. While we identified additional components of 
the Ha emission from a spiral density wave in the 
Be disk and from a second interaction region, the 
emission is not consistent with an accretion disk 



(McSwain et al. 2010). Nagae et al. (2006) also 



did not detect polarization from a second disk in 
LS I +61 303. High resolution X-ray spectroscopy 
should reveal an accretion disk from its thermal 
emission or a reflection iron line, but neither was 



detected in LS 5039 (Bosch- Ramon et al. 20071 



No high resolution grating spectra have been pub- 
lished for LS I +61 303, although an iron line is 



not required to fit low resolution spectra ( Harrison 
et al.[[2000[ [Chernyakova et al.[[2006[ ). 



Because the absence of an accretion disk is 
inconsistent with the microquasar model, |Dubus[ 
(2006) proposed that the extended radio emis- 



sion and observed X-ray/7-ray properties of LS 
I +61 303 and LS 5039 are better explained by 
an interaction shock between the stellar wind and 
a relativistic wind from a non-accreting pulsar. 



Dhawan et al. (2006) argue that VLB A images 
covering a full orbit of LS I +61 303 support the 
pulsar model rather than a precessing jet. Ha 
emission from LS I +61 303 also supports the 
wind shock model (McSwain et al. 20101. Such 



a model for LS I +61 303 was originally proposed 



by Maraschi & Treves ( 1981 ). The relativistic col 



liding wind mechanism successfully explains emis- 
sion from PSR B1259-63 ( [Cominsky et al.|1994[ ), 
the first radio pulsar discovered in a binary sys- 
tem with a main sequence companion (B2 Ve star; 
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Johnston et"al]|1992 1. 

Detecting a pulsar in LS I +61 303 or LS 5039 
would conclusively settle the debate over the high 
energy production mechanism. The pulsar timing 
solution would provide critical inputs for model- 
ing, such as the magnetic field strength and spin- 
down rate of the pulsar. It would also provide a 
tomographic probe of the circumstellar envelope, 
as demonstrated with PSR B1259— 63 ( ^Johnston] 
et al. 112005 1 . Finally, the pulsar timing would also 



improve the precision of the orbital parameters 
that have been measured via optical spectroscopy. 

Several studies have attempted to detect the 
putative pulsars, to no avail. LS 5039 lies within 
the strip of the Galactic plane surveyed by the 
Parkes Multibeam Pulsar Survey, which per- 
formed a low frequency search at 1374 MHz and 
flux density sensitivity 51374 = 0.14 mJy ( |Manch- 



ester et al. 2001). LS I +61 303 is too far north 



to be included in the Arecibo L-band Feed Array 
(ALFA) pulsar survey ([Cordes et al.||2006|, but 



it was included in the Green Bank Northern Sky 
Survey (Sayer et al. 1997). X-ray pulse searches 
in LS I +61 303 have been performed by[ 



et al. ( 1996 1 using RXTE, Harrison et al. 



with RXTE and ASCA, and by Rea et al. 



'aylor| 



(2000) 



(2010) 



with the Chandra ACIS-S camera, each with null 
results. 

To further investigate the possible presence of 
a pulsar in these systems, we performed a deep 
radio pulse search modeled after successful obser- 
vations of the radio pulsar PSR B1259— 63 near 



its periastron (Johnston et al. 2005). Our obser- 



vations were timed to occur near apastron, when 
the neutron star passes through a minimum in the 
surrounding stellar wind density. These observa- 
tions are described in §2. We model the stellar 
wind densities and optical depths, and discuss the 
null results from our pulsar search, in §3. 

2. Observations 

We observed LS I +61 303 and LS 5039 using 
the National Radio Astronomy Observatory's 100 
m Green Bank Telescope (GBT) and the Pulsar 
Spigot back-end (Kaplan et al.|2005) in mode 2. 



For LS I +61 303, we performed observations dur- 
ing 2008 May, September, and October using three 
receivers: S-band with 600 MHz of bandwidth 
(1650-2250 MHz), C-band with 800 MHz of band- 



width (4400-5200 MHz), and X-band, also with 
800 MHz of bandwidth (8500-9300 MHz). We 
observed LS 5039 using only the lower frequency 
S-band (1700-2400 MHz) and C-band receivers 
during 2008 September and November. Unfor- 
tunately we did not have sufficient time sched- 
uled to also include the X-band for LS 5039. The 
Spigot has 1024 available channels, but because 
of the smaller usable band at S-band, we only 
used 768 or 896 channels for those observations. 
The Spigot summed and synthesized the data into 
0.78125-MHz frequency channels every 81.92 fis. 
The times of mid exposure and the corresponding 
orbital phases are listed in Table [T] Illustrations 
of the orbital coverage of both systems are shown 
in Figures [T] and [2j 

The data were reduced using the PRESTO soft- 
ware package (Ransom 2001). We first searched 
the raw data for radio frequency interference in 
both the time and frequency domains and ap- 
plied interference masks to remove any interfer- 
ence. RFI was not a major issue with any of the 
searches except possibly for the S-band search of 
LS 5039, where approximately 20% of the data 
were recommended to be removed. For all the 
other observations the fraction was less than 2% of 
the data. As far as the fraction of the Fourier spec- 
trum is concerned, we removed less than 0.02% of 
the bins. LS I +61 303 and LS 5039 are expected 
to have interstellar dispersion measures (DM) of 
^ 100 cm-3 pc using the NE2001 model ( [Cordes 
& Lazio 2002 1 . We searched the observations 



by dedispersing the raw data into separate time 
series with DMs ranging from 0—2000 cm~'^ pc 
(to account for significant intersystem dispersion 
due to the dense stellar winds in both systems) 
and spaced by 0.5 cm~'^ pc. We then Fourier- 
transformed each time series and searched them 
by using Fourier-domain acceleration search tech- 
niques in order to maintain sensitivity to the pu- 
tative pulsars in their binary orbits. For the ac- 
celeration search, we search a range of frequency 
derivatives, which are assumed to be constant over 
the observation. This is an approximation to 
the true orbital acceleration of the binary system, 
which is sufficient for cases where the observation 
spans less than 10% of the orbital period, as is the 
case here. Since it is possible that pulsed emission 
can only been detected in brief visibility windows 
and the DM may change significantly during the 
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observations due to the clumped stellar winds, we 
also performed searches for single pulses using the 
standard tools in PRESTO. We searched for pulse 
frequencies up to 1 kHz in each case. 

A typical spectral in dex of pulsars is F^, ^ v~^-^ 
(Lorimer et al. 1995) or ~ v~^-^ (Maron et al 



2000), but the optical depth due to free- free ab- 
sorption by the stellar winds decreases as t^, ~ v~'^ 
( Lamers fc Cassinelli))1999p . Therefore the high 
frequencies used in our search improve the sensi- 
tivity of our observations by countering the effects 
of radio absorption and scattering by the compan- 
ions' winds. The flux density sensitivity limits are 
listed in Table [T] column 6. We did not detect any 
radio pulses from either LS I -|-61 303 or LS 5039. 

3. Discussion 

For the purpose of this discussion, we assume 
that the compact companions in both LS I -1-61 
303 and LS 5039 are 1.4 Mq neutron stars. Us- 
ing the mass functions, f{m), and their formal 
errors determined by Aragona et al. (20091, we 



can constrain the inclination of both systems. As- 
suming an optical companion with mass = 
12.5 Mq and radius = 6.1 Rq in LS I -^61 303 



(Casares et al. 2005a), the system inclination is 



i = 73°(+i7°^_9°). Likewise for LS 5039, we as- 
sume an optical companion of 28.5 Mq and 9.7 Rq 



based on the spectral type ( Martins et al. 2005 1 , 
and i = 71°(+19°,-8°). 

The wind velocity, v{r), of massive stars can be 
reasonably well estimated by the relation 



v{r) 



1 - 



R. 



(1) 



(Lamers & Cassinelli 1999). For an 0-type star 
with weak winds, as is LS 5039, the coefficient 



(3=1 (Puis et al. 1996). Stochastic clumping is 



ubiquitous among the winds of O stars, and we 
quantify these clumping effects using the method 



described by |Puls et al.| ([2008j) . The mean stellar 
wind density is 



(p(0> 



M 



(2) 



47r v{r) 

We assume that the space between clumps is a 
vacuum, and the density within the clumps, pc, is 
enhanced by a volume filling factor /vol such that 

(P> = /volPc- (3) 



The absorbing path length is reduced accordingly 
since no absorption occurs in the inter-clump re- 
gion. For LS 5039, we assume a spherically sym- 
metric wind. LS I -1-61 303 contains a rapidly 
rotating Be star, and below we discuss how the 
structure of the winds may be modified by the 
fast rotation. 

LS 5039 has a terminal wind velocity, Voo, 
of 2440 km s~^. Its mass loss rate is M « 
10~^ Mq yr~^, although this rate is variable by a 
factor of 2 ( [McSwain et~aI1|2004[ ). LS 5039 is not 
known to be a rapidly rotating star, so no signifi- 
cant deviation from spherically symmetric winds is 
expected. Assuming a fully ionized hydrogen wind 
with a temperature of 2/3 the effective tempera- 
ture of the O star (Toff = 37500 K; [McSwain et i 
2004|, the optical depth, t^, ~ due to free-free 
absorption of the resulting wind density is given 
by Lamers & Cassinelli (1999). Using this simple. 



spherically symmetric wind model for LS 5039, we 
integrated along the line of sight to the putative 
neutron star's position in the inclined orbit to de- 
termine the free-free optical depth as a function 
of orbital phase. Since the wind is a highly ion- 
ized plasma, other scattering and/or absorption 
mechanisms may also be significant: pulses may be 
smeared due to the variable column density or tur- 
bulence in the clumpy wind, nonlinear optical pro- 
cesses may scatter photons due to the high photon 
occupation number in the vicinity of a pulsar, and 
absorption at the cyclotron resonance frequency 
and its harmonics may be significant (jThompson 
et al.|1994 ) . We neglect these additional processes 



here. 

Our resulting for LS 5039 are plotted in Fig- 
ure [s] for /vol = 1 (corresponding to an unclumped 
wind). Since the free- free opacity has a depen- 
dence, the effect of wind clumping is to enhance 
Ti, by a factor of l//voi despite the reduced path 
length of absorbers. Evolved, Galactic O stars 
have 0.1 < /vol < 0.5 (Puis et al.|[2008| and ref- 



erences therein), although the less evolved, main 
sequence O star in LS 5039 may have less strongly 
clumped winds. Since LS 5039's mass loss rate is 
variable, the predicted optical depth may also vary 
by a factor of ^ 2 at the time of our observations. 
The large ^ 1 are prohibitive for detecting any 
possible pulsar in the LS 5039 system at the S- 
band and C-band frequencies used in our search. 
However, we have neglected the effects of a con- 
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ical wind shock region, which hkely shields the 
putative pulsar in LS 5039 from the dense stel- 
lar wind. If the O star's winds are colliding with 
a relativistic pulsar wind, then the wind collision 
zone would result in a double shock structure sep- 



with latitude (Maeder 2009), leading to a slower 



arated by a contact discontinuity (CD; Bogovalov 



et al. 20081. The structure has a conical shape 
that envelopes the compact object or the hot star, 
depending on the ratio of momentum flux of the 
winds, 1] = {E/c)/{Mv). This ratio uniquely de- 
termines the opening angle and location of the 
shock structure, and it will define the topography 
of a wind cavity bounded by the CD. At the time 
of inferior conjunction, the absorption along the 
neutron star's line of sight may be far less than 
predicted for a uniform stellar wind. High reso- 
lution, high signal-to-noise, phase-resolved obser- 
vations of the ultraviolet wind resonance lines are 
necessary to determine how large a cavity is evac- 
uated within the stellar winds given the fast orbit 
of LS 5039. 

The wind properties of LS I +61 303 are un- 
certain, but we estimated its M k, 10~® Mq yr~^ 



using the method of Lamers & Cassinelli (1999) 
for B-type stars. The optical star in this sys- 
tem is a rapidly rotating Be star, probably ro- 
tating near the critical velocity at which gravity 
balances the centrifugal force at the equator. For 
a typical Be star, Kot/K:rit ~ 0.8 (McSwain et 
ar]|2008 ). The fast rotation leads to several inter- 



esting effects at the stellar surface that affect the 
mass loss. Due to the large centrifugal force at the 
equator, the Roche model predicts that the equa- 
torial radius i?oq is about 1.3 times larger than 



the polar radius i?p ( Maeder||2009 ). According to 
the von Zeipel theorem, the smaller effective sur- 
face gravity at the equator implies a cooler equa- 
torial temperature, Teq_eff, and decreased equato- 
rial flux ( |von Zeipel||1924| . We modeled LS I +61 
303 using a polar temperature of Tp off — 30, 000 
K and T^q^cS = 21,500 K, with a mean surface 
temperature Tcff — 25, 100 K. This Tgff is consis- 



tent with measurements by Howarth ( 1983 ) and 
references therein. The cooler Toq^off leads to an 
increased opacity, enhancing the equatorial mass 



loss rate despite the lower radiative flux (Maeder 
2009 ). In Figure |4j we plot the resulting mass flux 



as a function of stellar latitude, normalized to a 
total M = 10^^ Mq yr^^. The wind momentum 
is also affected by the variation of radiative flux 



wind velocity at the equator (also shown in Figure 
[3]). Combining these effects, we find a significantly 
denser wind at the equator, so our model is a rea- 
sonable representation of the dense circumstellar 
disk of the Be star. Finally, we assume that the 
equatorial disk is aligned with the binary orbit, 
we neglect wind clumping effects, and we integrate 
through the resulting wind distribution along the 
line of sight to the putative pulsar to calculate 
during the orbit. 

We plot the resulting as a function of or- 
bital phase for LS I +61 303's asymmetric wind 
in Figure [5] Only the X-band provides an op- 
portunity for an optically-thin view of the pulsar 
during 0.4 < (/)(TG) < 0.8, assuming the system 
inclination is i = 73°. Due to telescope schedul- 
ing constraints, our GBT X-band observation was 
performed at (/'(TG) = 0.97. While our model sug- 
gests that Tjf < 1 for the lower limit of i, it is not 
an optimal phase for detecting pulses. If the incli- 
nation is high, the pulsar will be viewed through 
the densest portion of the wind and would likely 
be highly absorbed at all orbital phases. 

While the Be star's disk is truncated to a max- 
imum radius of 34-37 Rq (Grundstrom et al. 



2007), we predict that the pulsar can be eclipsed 
by the disk if the orbit is viewed edge-on, lead- 
ing to very high near 0(TG) =0.1. However, 
in this post-supernova system, there is no require- 
ment that the orbital plane remain aligned with 
the circumstellar disk. A kick velocity from the su- 
pernova would likely have disrupted the system's 
angular momentum sufficiently to misalign the or- 
bit. Indeed, the Be star disk in PSR B1259-63 is 



highly inclined relative to the orbital plane ( John 



ston et al. 


1999 


Bogomazov 


2005 



Therefore 

a pulsar in LS I +61 303 may not be strongly 
eclipsed at any phase. Further observations to 
constrain both the stellar inclination as well as 
the orbital inclination are necessary to improve 
our model of the stellar wind density of LS I +61 
303. 

Based on our model, we expect that the free- 
free absorption is prohibitively large at the S-band 
and C-band frequencies to detect a pulsar at any 
orbital phase of either LS I +61 303 or LS 5039. 
The high frequency X-band offers the best chance 
of countering the effects of stellar wind absorption. 
Although we did not detect pulsars in this search. 
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we conclude that detecting a pulsar in LS 5039, 
using higher frequencies, may be possible if the 
CD due to a colliding wind shock region effectively 
shields the pulsar from the dense stellar wind. De- 
tecting a pulsar in LS I +61 303 should be possible 
using the X-band during a more optimal orbital 
phase as long as the pulsar beam crosses our line 
of sight. 

In order to power these 7-ray sources, any asso- 
ciated pulsar should have a high spin-down power 
> 10^^ erg s~^. Such a pulsar likely has a rota- 
tion period ~ 0.05-0.4 s, with an average beaming 



fraction of 20-50% (Tauris & Manchester 19981 



There are only four known 7-ray binaries with well 
established MeV-TeV emission: LS I -f61 303, 
LS 5039, PSR B1259-63, and HESS J0632-F05 7 
([Aharonian et al. 2007 Bongiorno et al. 20111. 
PSR B1259— 63 remains the only system with a de- 
tected radio pulsar, but with further observations 
it may be possible to expand this elite group. 

Future X-ray pulsation searches of these sources 
are also encouraged. We note that the HMXB 4U 
0115-1-63 (= V635 Cas) has a very similar spectral 
type and orbital period compared to LS I -1-61 303 
(B0.2 Ve, 24.3 d in this case; ( |Liu et al.|r2006[ )) 
and their stellar wind properties should be very 
similar. That system is known to include an ac- 
creting pulsar (e.g. Nakajima et al. 2010). The 



thermal emission of the hot star winds, usually 
with 0.2 < kT < 0.6 keV but sometimes hotter 



(Naze et al. 2011) will likely prohibit the detec- 



tion of thermal surface emission from a neutron 
star. Harder, non-thermal X-rays (above 2 keV) 
should be more effective for pulsation searches in 
these 7-ray binaries. 

We are grateful to Carl Bignell and the staff 
at NRAO for their help scheduling and perform- 
ing these observations. We also thank Guil- 
laume Dubus and Anna Szostek for helpful dis- 
cussions about the stellar winds and possible con- 
tact discontinuities in these systems. The Na- 
tional Radio Astronomy Observatory is a facil- 
ity of the National Science Foundation operated 
under cooperative agreement by Associated Uni- 
versities, Inc. This work is supported by NASA 
DPR numbers NNX08AV70G, NNG08E1671, and 
NNX09AT67G. MVM is grateful for an institu- 
tional grant from Lehigh University. 
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Fig. 1. — Orbital phase coverage of our GBT ob- 
servations of LS I +61 303 (gray) compared to the 
orbital geonietry (Aragona et al.|[200 9). The rela- 
tive orbits (r/a) of the 12.5 Mq optical star (dot- 
ted line; Casares et al. 2005a[ ) and a presumed 
1.4 Mq neutron star companion (solid line) are 
shown. The center of mass is indicated with a 
cross, and the thin solid line is the orbital major 
axis. Points around the orbit indicate steps of 0.1 
in orbital phase, and the phases of periastron and 
apastron are labeled. 
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Fig. 2. — Orbital phase coverage of our GBT 
observations of LS 5039 (gray) compared to the 
orbital geometry, in the same format as Fig. [l] 
Here, the optical star (28.5 Mq) orbits a presumed 
1.4 Mq neutron star. 
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Fig. 3. — The calculated optical depth, r^, of LS 
5039 as a function of orbital phase. Errors in 
due to the possible range in i are shown in gray. 
The orbital phases covered by our observations are 
marked for each receiver. The calculation assumes 
a spherically symmetric wind but neglects the ef- 
fects of an interacting wind shock region. 



Fig. 4. — The mass flux (top) and terminal wind 
velocity Voo (bottom) for LS I +61 303, assuming 
a rotationally distorted stellar wind. The cooler 
7oq,off leads to decreased equatorial radiative flux 
but increased opacity, enhancing the equatorial 
mass loss rate. The mass flux is normalized to 
a total M = 10"^ Mq yr'^ The lower radiation 
fltix decreases the wind momentum and thus Voo 
at the equator. 
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Table 1 
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MJD of 


Orbital 


Exposure Time 




'S'lim 


btar 


Mid Exposure 


Phase 


(hr) 


Receiver 




LS I +61 303 


54599.631 


0.96 


3.0 


C-band 


4.1 




54599.777 


0.97 


3.0 


X-band 


7.5 




54599.882 


0.97 


1.2 


S-band 


9.4 




54729.006 


0.85 


0.5 


S-band 


14.5 




54760.010 


0.02 


1.5 


C-band 


5.8 


LS 5039 
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Fig. 5. — The calculated optical depth, Ti,, of LS I 
+61 303 as a function of orbital phase, using the 
ephemeris of Taylor & Gregory (1982 1. Errors in 
due to the possible range in i are shown in gray 
for the X band only; the C and S band errors are 
comparable. The orbital phases covered by our 
observations are marked for each receiver. The 
calculation assumes a rotationally distorted star 
whose equatorial plane is aligned with the orbital 
plane, and it neglects the effects of an interacting 
wind shock region. 
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